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Abstract — Multi-antenna precoding effectively mitigates the 
interference in wireless networks. However, the precoding ef- 
ficiency can be significantly degraded by the overhead due to 
the required feedback of channel state information (CSI). This 
paper addresses such an issue by proposing a systematic method 
of designing precoders for the two-user multiple-input-multiple- 
output (MIMO) interference channels based on finite-rate CSI 
feedback from receivers to their interferers, called cooperative 
feedback. Specifically, each precoder is decomposed into inner 
and outer precoders for nulling interference and improving the 
data link array gain, respectively. The inner precoders are further 
designed to suppress residual interference resulting from finite- 
rate cooperative feedback. To regulate residual interference due 
to precoder quantization, additional scalar cooperative feedback 
signals are designed to control transmitters' power using different 
criteria including applying interference margins, maximizing sum 
throughput, and minimizing outage probability. Simulation shows 
that such additional feedback effectively alleviates performance 
degradation due to quantized precoder feedback. 

I. INTRODUCTION 

In multi-antenna wireless networks, precoding can effec- 
tively mitigate interference between coexisting links. This 
paper presents a new approach of efficiently implementing pre- 
coding in the two-user multiple-input-multiple-output (MIMO) 
interference channels by exchanging finite-rate channel state 
information (CSI). Specifically, precoders are designed to 
suppress interference to the interfered receivers based on 
their quantized CSI feedback, and the residual interference 
is regulated by additional feedback of power control signals. 

Recently, progresses have been made in analyzing the 
capacity of multi-antenna interference channels. In particular, 
interference alignment techniques have been proposed for 
achieving the channel capacity for high signal-to-noise ratios 
(SNRs) [1]. Such techniques, however, have limited practi- 
cality due to their complexity, requirement of perfect global 
CSI and their sub-optimality for finite SNRs. This prompts 
the development of linear precoding algorithms for practical 
decentralized wireless networks [2]-[5]. For time-division 
multiplexing (TDD) multiple-input- single-output (MISO) in- 
terference channels, it is proposed in [2], [5] that forward- 
link beamformers can be adapted distributively based on 
reverse-link signal-to-interference-pluse-noise ratios (SINRs). 
Targeting the two-user MIMO interference channels, linear 
transceivers are designed in [3] under the constraint of one 
data stream per user and using different criteria including 
zero-forcing and minimum mean square error. In [4], the 
rate region for MISO interference channels is analyzed based 
on the cognitive radio principle, yielding a message passing 



algorithm for enabling distributive beamforming. The above 
prior work does not address the issue of finite-rate feedback 
though it is widely used in the practice to enable precoding. 
Neglecting feedback CSI errors in precoder designs leads to 
over optimistic network performance. 

For MIMO precoding systems, the substantiality of CSI 
feedback overhead has motivated extensive research on CSI 
quantization algorithms, forming a field called limited feed- 
back [6]. Recent limited feedback research has focused on 
MIMO downlink systems, where multiuser CSI feedback 
supports space division multiple access (SDMA) [7]. It has 
been found that the number of feedback bits per user has to 
increase with the transmit SNR so as to bound the throughput 
loss caused by feedback quantization [8]. Furthermore, such 
a loss can be reduced by exploiting multiuser diversity [9], 
[10]. Designing limited feedback algorithms for interference 
channels is more challenging due to the decentralized network 
architecture and the growth of feedback CSI. Cooperative feed- 
back algorithms are proposed in [11] for a two-user cognitive 
radio network, where the secondary transmitter adjusts its 
beamformer to suppress interference to the primary receiver 
that cooperates by feedback to the secondary transmitter. This 
design is tailored for a MISO cognitive radio network and thus 
unsuitable for general MIMO interference channels, which 
motivates the current work. 

We consider two coexisting MIMO links where all nodes 
employ equal numbers of antennas and linear precoding is 
enabled by quantized cooperative feedback. Channels are 
assumed to have i.i.d. Rayleigh fading. A systematic method is 
proposed for jointly designing linear precoders and equalizers 
under an orthogonality constraint, which decouples the links 
in the case of perfect feedback. To be specific, precoders and 
equalizers are decomposed into inner and outer components, 
where the former are designed to suppress residual interference 
caused by feedback errors and the latter to enhance link array 
gain. Second, additional scalar cooperative feedback, called 
interference power control (IPC) feedback, is proposed for 
controlling transmitters' power so as to regulate the residual 
interference. Specifically, the IPC feedback algorithms are 
designed using different criteria including fixed interference 
margin, maximum sum throughput, and minimum outage 
probability. 

Notation: The superscript f represents matrix Hermitian 
transpose. The operator [X]& gives the fcth column of a matrix 
X. Let -<, >: and y represent element- wise inequalities 
between two real vectors. 
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Fig. 1. MIMO interference channels with local and cooperative feedback. 

II. System Model 

We consider two interfering wireless links as illustrated 
in Fig. 1, where the two pairs of transceivers are denoted 
as (Ti,Ri) and (T2, R2). Each transmitter/receiver has L 
antennas employed for suppressing the interference as well 
as supporting spatial multiplexing. These functions require 
CSI feedback from receivers to their interferers and intended 
transmitters, called cooperative feedback and local feedback, 
respectively. We assume perfect CSI estimation and local 
feedback, allowing the current design to focus on suppressing 
interference caused by cooperative feedback quantization. All 
channels are assumed to follow independent blocking fading. 
The channel coefficients are samples of i.i.d. £/V(0, 1) random 
variables. Let H mn denote a L x L i.i.d. £A/"(0, 1) matrix 
representing fading of the channel from T n to R m . Then the 
interference channels are modeled as {z/H mn } and the data 
channels as {H mm } where m,n G {1,2} and m ^ n. The 
factor v < 1 quantifies the difference in transmission distance 
between the data and interference links. 

Each link supports M < L spatial data streams by linear 
precoding and equalization. To regulate residual interference 
caused by precoder feedback errors, the total transmission 
power of each transmitter is constrained by cooperative IPC 
feedback. For simplicity, the scalar IPC feedback is assumed 
perfect since it requires much less overhead than the precoder 
feedback. Each transmitter uses identical transmission power 
for all spatial streams, represented by P n for T n with n = 1, 2 
and the maximum P max . Assume that all additive white noise 
samples are i.i.d. CA/"(0, 1) random variables. Let G m and 
F m denote the linear equalizer used by R m and the linear 
precoder applied at T m , respectively. Thus the receive signal- 
to-interference-plus-noise ratio (SINR) at R m for the ^th 
stream can be written as 

SINRfg := P ™\[G m ]\K [¥ m ] e \ 2 m + n _ (J) 
l + P„H|[G ro ]]H TOn F„||2 

Two performance metrics, ergodic throughput and outage 
probability, are considered. The total ergodic throughput of 



both links, called sum throughput, is defined as 

2 M 

C:=££E[log 2 (l + SINRM) 



(2) 



where SINR^J is given in (1). Next, consider the scenario 
where the coding rates for all data streams are fixed at 
log 2 (l + 0) where 6 is the receive SINR threshold for correct 
decoding. We define an outage event as one that the SINR of 
at least one data stream is smaller than 6. It follows that the 
outage probability is given by 



Pont := Pr 



min min 

n=l,2 1<£<M 



(3) 



III. Precoding with Limited Feedback 

A. Precoder Design 

A pair of precoder and equalizer (G m , F n ) with m ^ n are 
jointly designed under the following orthogonality constraint 



m, n G {1, 2}, m ^ n. 



(4) 



The constraint aims at decoupling the links and requires that 
L > 2M. A key step of the proposed design is to decompose 
the precoder F n into an inner precoder and an outer 
precoder F° . Specifically, F n = F^F° where F^ and F° are 
L x M and M x M matrices, respectively, where the size of 
F^ is minimized to reduce feedback overhead. Similarly, we 
decompose the equalizer G rn as G rn = G' m G°, where G' m 
is a L x N inner equalizer and G^ aiVxM outer equalizer, 
where N is a design parameter under the constraints N > M 
and N < L — M. The inner precoder/equalizer pair (G^, F^) 
is designed to enforce the constraint in (4) while the outer 
pair (G^,F°) enhances the link array gain as discussed in 
the sequel. It follows that 



(G^) f H mn F^ = 0, m^n. 



(5) 



Under this constraint, (G^, F^) are designed by decomposing 



FL 



using the singular value decomposition (SVD) as 
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where the unitary matrices V mn and U mn consist of the left 
and right singular vectors of H mn as columns, respectively, 
and jAl5 n j denote the eigenvalues of H mn Hj nn following 
the descending order. Let A and B be two subsets of the 
indices {1, 2, • • • , L} with \A\ = N, \B\ = M, and AnB = 0. 
The constraint in (5) can be satisfied by choosing 

GL = {[V mn ]fc I k G A} and F^ = {[U mn ] fc | k G B}. (7) 

Given (G^FjJ, the outer pair (G^,F^) are jointly 
designed based on the SVD of the N x M effective channels 
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(8) 



Note that the elements of H^ m are i.i.d. CA/*(0, 1) random 
variables and their distributions are independent of (G^, F^J 
since H mm is isotropic. To transmit data through the eigen- 
modes of H mm , G^ and are chosen as 

G°, = {[V mm ]i, [V mm ] 2 , • • • , [V mm ]M} and F^ = U mm . 

With perfect CSI feedback, the above precoder and equalizer 
joint design converts each data link into M parallel spatial 
channels which are free of interference. 

Note that increasing TV enhances the array gain of both 
links. Specifically, the expectations of the SNRs increase with 
TV. Thus, TV should take its maximum (L — M). However, 
maximizing TV need not be optimal for the link performance 
in the case of quantized feedback as discussed in the sequel. 

B. Quantized Precoder Feedback 

In this section, we choose the index sets A and B in (7) with 
the objective of suppressing the residual interference caused 
by precoder feedback errors. 

Recall that the precoding at T n is enabled by quantized 
cooperative feedback of F^ from R m with m ^ n. Let F^ 
denote the quantized version of F^ and define the resultant 
quantization error e n as 



(Fk) f FU| 2 F 
M 



1,2 



(9) 



where < e n < 1. The error e n is zero in the case of 
perfect cooperative feedback, namely F^ = F^. A nonzero 
error results in violation of the orthogonality constraint in (5) 



G^H mn F^ ^ 0, m^n. 



(10) 



The resultant residual interference from T n to the data 
stream of R m has the power 



^mn :— ^n^|| [G m ]J(G m )''"H mn F n F n || , 



m^n. (11) 



Next, we choose A and B in (7) by minimizing an upper 
bound on the residual interference power as follows. Based on 
(7), (6) can be rewritten as 



where T, r 



[B m G^J 5] mn [F^ C n (12) 

i diag fvS, V5i> • • • > 4^ 



with II being an arbitrary permutation matrix that rearranges 
the order of the singular values along the matrix diagonal. The 
columns of the matrices B m and C n comprise the (L — TV) 
left and (L — M) right singular vectors of H mn , respectively, 
which are determined by II. Let the set V mn contain the last 
TV elements along the diagonal of E mn . 



Lemma 1. The interference power I^n in (11) can be upper 
bounded as 

I$ n < MvP n e n max a 2 , m 7^ n. (13) 

aeT> mn 

Readers can refer to the full paper [12] for the proofs of 
the above lemma as well as analytical results in the sequel. 
Minimizing the upper bound in (13) gives that V mn consists 
of the TV smallest singular values of H mn . Equivalently, II is 
an identity matrix and thus G^ and F^ are given as 



Gm = {[V mn ]L-JV+l, [V mn ]L-JV+l, • • • , [V mn ]z,} 



(14) 



Fn — {[U mn ]i, [U mn ]2, • • • , [U mn ]M} • 

Then (13) can be simplified as 

4i < MvP n \%- N+1 h n , Vl<KM,m/n. (15) 

Note that the above upper bound on linn is independent of 
the stream index L On one hand, the upper bound reduces 
with decreasing TV. On the other hand, as mentioned earlier, 
larger TV increases link array gain. These opposite effects of TV 
on link performance make it an important parameter for pre- 
coder optimization. Finding the optimal TV is mathematically 
intractable but a numerical search is straightforward. 

IV. Interference Power Control Feedback 

A. Fixed Interference Margin 

The receiver R m sends the IPC signal, denoted as r] n , to the 
interferer T n for controlling its transmission power as 



P n = min(^ n ,P max ), 



1,2. 



(16) 



The scalar r] n is designed to prevent the per-stream interfer- 
ence power at R m from exceeding a fixed margin r with r > 0, 
namely I$ n < r for all < i < M. A sufficient condition 
for satisfying such constraints is to upper bound the right hand 
side of (15) by r. It follows that 

~ .[L-iV+l] > m + U ' W 
Mv\\nn 

Given r, a lower bound Aim on the sum throughput C, called 
the achievable throughput, is obtained from (2) as 



M 



m=l x 



min(r/ m ,P m ax )A 

1+7 



(18) 



where 7? m is given (17). 

It is infeasible to derive the optimal value of r for either 
maximizing A\y\ in (18) or minimizing P out in (3). However, 
for P max being either large or small, simple insight into 
choosing r can be derived as follows. The residual interference 
power decreases continuously with reducing P max . Intuitively, 
r should be kept small for small P m ax- For large P ma x, the 
choice of r is less intuitive since large r lifts the constraints 
on the transmission power but causes stronger interference and 
vice versa. We show below that large r is preferred for large 
P max . Let Afc denote the eigenvalue of the Wishart matrix 
HHt with H being an i.i.d. TV x M CA/*(0, 1) matrix. Define 
Afc similarly but with H being aLxL matrix. 



the achievable throughput is 



(1 + r)M^A L _Ar+iei 



+ o(l). 



Lemma 2. For large P ma: 

Ml" / 

A M =2^E log 2 11 + 

t=i L V 

It can be observed from the above result that the first order 
term of Aim attains its maximum for r — >> oo. However, this 
term is finite even for asymptotically large P ma x and r, which 
is the inherent effect of residual interference. 



Lemma 3. For large P m 
bounded as 

Pout<2Pr 



the outage probability is upper 



<e\ +o(i). 



x l + r M^A L _Ar+iei 
Similar remarks on Lemma 2 apply to Lemma 3. 

B. Sum Throughput Criterion 

In this section, an iterative IPC algorithm is designed for 
increasing the sum throughput C in (2). Since C is a non- 
convex function of transmission power, directly maximizing 
C does not yield a simple IPC algorithm. Thus, we resort to 
maximizing a lower bound A$j (achievable throughput) on C 
instead, obtained from (2) and (15) as Ast = IE [A] with 



M 
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log 2 1 + 



p \ * 



l + P 1 Mv\ [ 2 L - N+1] e 1 , 



(19) 



Thus, the optimal transmission power pair is given as 

(PZ,P1)= max A{P 1 ,P 2 ). 

Pi,P2e[o,p m 



(20) 



The objective function A remains non-convex and its max- 
imum has no known closed-form. However, inspired by the 
message passing algorithm in [4], a sub-optimal search for 
(Pj*, P 2 ) can be derived using the fact that 

dA(p^p^) 



OPrr 



V ra = 1,2. 



To this end, the slopes of A are obtained using (19) as 

dA(P u P 2 ) 



OPrr 



Pn 



Pm 



(21) 



where 



log 2 



EM 
I 1 i 



)+ 
A r< 



^ m := log 2 e 



M 



+ Mv\tn N+1] e n P n - 

i [L-N+l] £ 



/ ^rara r m 



A/fry \ L ^ _iV ^ ±J ^ P J- \M P 
iK/ is /\nm ^rn r rn \ / ^nn r n 



log 2 eM 2 ^A^ N+1] e m 
Pm I_ 1 + Mu\^ L - N+1 ^ P ' 

1 -|- iVi VA nrn e m .r m 

Note that based on available CSI, /i m has to be computed at 
R m and (ip m ,Pm) at P n with n ^ m. Therefore, based on 
(21), an iterative IPC feedback algorithm can be designed to 
have the following procedure. 



Algorithm 1: 

1) The transmitters Ti and J 2 arbitrarily select the initial 
values for Pi and P2, respectively. 

2) The transmitters broadcast their choices of transmission 
power to the receivers. 

3) Given (Pi,P 2 ), the receiver Ri computes (^1,^2^2) 
and sends /ii and (ip2 — P2) to Ti and T 2 , respectively. 
Likewise, R2 computes (p2-,^i- ) pi) an d feeds back /12 
and (^1 — pi) to T 2 and Ti, respectively. 

4) The transmitters Ti and T 2 update Pi and P 2 , respec- 
tively, using (21) and the following equation 



Pm(k+1) = min < 



where k is the iteration index and A7 a step size. 
5) Repeat Steps 2) — 4) till the maximum number of 
iterations is performed or the changes on (Pi,P 2 ) are 
sufficiently small. 

C Outage Probability Criterion 

As the problem of minimizing P out in (3) by power control 
is analytically intractable, the IPC algorithm is designed by 

minimizing an upper bound on P out . Using (15), the SINR in 

M 

(1) is lower bounded by SINR m where 



SINR 

Therefore, 
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m /x mm 



1 + P n Mv\tn N+1] e r , 



m ^ n. 



Pout < 



Pr 



min SINR 

71=1,2 
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(22) 



(23) 



Minimizing the above upper bound on P out is similar in 
the mathematical structure to the classic problem of optimal 
power control for single- antenna interference channels [13]. 
The optimal transmission power for minimizing the right hand 
side of (23) solves the following optimization problem 



arg mm 

Pi,p 2 e[o,p max 



S.t. min SINR 

m=l,2 



[M] 



(Pi,P 2 )>0 



(24) 



where the first minimization implies that (Pj*, P 2 ) ^ (Pi, P2) 
for all (Pi,P 2 ) that satisfies the constraint in (24) as well as 
(0,0) r< (Pi,P2) ^ (Pmax, Pmax), called feasible power pairs. 
Using (22) and (24), the constraint in (24) can be written as 

Pfn 

> a m + m ^ n (25) 

MuXt- N+1] e n 9 



where a Tl 



7[my and 



AM] 



power (A, A) that satisfies the constraints in (25) is 

' <2 m + h nin d n 



The minimum 



bmnbnm 



m 



1,2. 



(26) 



This expression gives optimal power control as stated below. 

Proposition 1. // (Pj\ P 2 * ) i/i (24) gxwto, (Pf, P 2 *) = (A, A) 
wzY/z (Pi,P2) given in (26). 
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Fig. 2. Comparison of achievable sum throughput between different IPC 
feedback algorithms for the coupling factor v = 0.2. 

Based on Proposition 1, the corresponding IPC feedback 
procedure is obtained as follows. 
Algorithm 2: 

1) The receiver Ri computes bi2 and transmits bi2 to T 2 . 
Similarly, R2 computes 621 and feeds back 621 to Ti. 

2) The receiver Ri computes a\ and communicates a\ to 
Ti. Likewise, R2 computes CL2 and feeds back CL2 to T 2 . 

3) The transmitters Ti and T2 compute P\ and P2, respec- 
tively and set their transmission power equal to (Pi, P2) 
if they are feasible. Otherwise, arbitrary transmission 
power is used. 

V. Simulation Results 

In the simulation, codebooks for quantizing F\ and F2 
are randomly generated and have equal sizes. The simulation 
parameters are set as L = 6, M = 2, N = 3, and B = 6. 
The interference margin is either fixed at r = 2 or increased 
as r = 0.4P max . 

Fig. 2 compares the achievable throughput of different 
IPC feedback algorithms. Significant coupling {y = 0.2) 
between links is observed to decrease achievable throughput 
dramatically with respect to perfect CSI feedback. For large 
^Pmax, the IPC feedback Algorithm 1 designed for maximizing 
the achievable throughput is observed to provide substantial 
throughput gain over those based on interference margins. 
Furthermore, increasing r with growing P max gives higher 
throughput than fixed r, which is consistent with Lemma 2. 

Fig. 3 compares the outage probabilities and average trans- 
mit SNRs of different IPC feedback algorithms. With respect 
to the IPC feedback using increasing r or with perfect feed- 
back, Algorithm 2 dramatically decreases average transmission 
power. Moreover, Algorithm 2 yields lower outage probability 
than the two algorithms using t. Finally, fixing t causes Pout 
to saturate as P ma x increases. 
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